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A theory of exciton condensed phase creation in two-dimensional system is presented. The theory 
is applied to explain the appearance of the periodical fragmentation which was observed last years 
in luminescence from the ring around laser spot in crystal with double quantum wells. 
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An electrons and holes of indirect excitons in semiconductor coupled quantum wells are localized 
in different wells and as the result the excitons have long lifetime. Such system present interest 
for study of exciton condensation phenomena jT|. Recently interesting experiments were fulfilled for 
indirect excitons in AlGaAs and InGaAs crystals with coupled quantum wells. In works [21 El S0 El El 
the authors have observed a ring structure in emission of indirect excitons outside of laser spot. The 
distance between ring and center of the excitation spot growths with increasing the pumping and 
can reach several hundred microns that is much larger than the exciton diffusion length. At low 
temperature (T < 2K) the external ring is fragmented into periodical structure over macroscopic 
lengths. The mechanism of the luminescence ring formation was proposed in [HI [Z] - The origin 
of fragmentation along the ring is unclear. In present work the appearance of this fragmentation is 
explained by creation of islands of exciton condensed phases along the ring. In [313 it was shown that 
in the presence of an attractive interaction between excitons a high density exciton system is unstable 
with respect to the periodical structure appearance. To obtain the spacial dependence of the exciton 
generation rate we shall use the suggestion of paper [HI Ej that in absence of the light irradiation the 
quantum well contains electrons and that holes created by light are captured by quantum well more 
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effectively than electrons. As the result the charge separation takes place from positive charge value 
in the center of laser spot till negative one far from the spot. The emission of light occurs from a ring 
where product of the electron and hole densities has maximum. For determination of the density 
distribution we have used the system of equations similar to that investigated in [DJ [7] but we added 
the equation for the exciton density. Lets n e , n h ,n ex designate the well densities of electrons, holes 
and excitons, respectively. These values satisfy to the following kinetic equations 

°"' D e A 2 n e + K e (r)-Wn e n h -^—^, (1) 



at t, 



e 



^ = D h A 2 n h + K h (r) - Wn e n h - ^, (2) 
at r h 

^ = D ex A 2 n e + G — — , (3) 
dt r ex 

where D e , Dh and D ex are the diffusion coefficients, r e , and r ex are the lifetime for electrons, holes 
and excitons, respectively, W is the electron-hole recombination rate, K e (Kh) is the electron (hole) 
creation rate in the well, G is the exciton production rate {G = qWn e nh, q < 1). 

Fig.l shows the result of calculations of electron, hole and exciton densities for the following 
values of parameters: D e = 200cm 2 /s, = 50cm 2 / s,D ex = 10cm 2 /s, W = 150cm 2 /s , q = 0.9, 
r e = Th = 10~ 5 s, r ex = 10 _7 s, the spacial distribution of the pumping was approximated by the 
Gaussian curve with width 60/im and Kh = 2i^ e . The qualitative picture of behavior of 
function of r is similar to results obtained in [7j. The position of the exciton density maximum 
determines the ring luminescence in experiments. We shall use the obtained value of exciton density 
to study the fragmentation of the exciton condensed phase. It is suggested that the condensation 
occurs due to the exciton-exciton interaction and it is not the Bose-Einstain condensation in the 
wave vector space with k=0. We have estimated the Van der Waals attraction using the well known 
formula of quantum mechanics and showed that Van der Waals interaction between two indirect 
excitons in coupled quantum wells exceeds the dipole-dipole interactions in the range of distances 
between excitons of order of several exciton radiuses (See Appendix). The attractive interaction may 
cause the existence of the exciton condensed phase. Some models of phase transitions in a system of 
indirect excitons were studied in [ID] . Afterwards we shall suggest that condensed phase exists and it 
is described by some parameters which will be determined later. Due to a finite value of the exciton 
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Figure V. The dependence of electron, hole and exciton densities on distance from center 

lifetime the sizes of exciton condensed phases are restricted. As the result in the two-dimensional 
case the condensed phase must exist as a system of islands similar to electron-hole drops in bulk 
semiconductors. In studied system the islands should be localized on the ring where the exciton 
density has maximum. 

Under following consideration we shall apply the theory of creation of the exciton condensed 
phase islands in two-dimensional case [11 J to studied system with nonuniform pumping. Firstly, we 
shall study the system of condensed phase islands periodically situated along the ring at maximum 
of the exciton density with r = tq (Fig. 2). (p m = 27nn/N is the angle of m-th island on the ring. 
Afterwards we shall show that deviations (fluctuations) from periodicity are small. Let us consider 
the formation of some island, for example,the island with m = 0. We introduce the distribution 
function /„, which determines a probability of the island with m = to have n excitons. The size of 
disks is determined by four processes: creation of excitons by pumping, capture of the excitons from 
environment, escape of excitons from the disk, and the exciton decay. The kinetic equation for the 
distribution function have the following form 

-gj- = -Jn+l + Jm (4) 

where j n is the probability current 

j n = 2irR n _ 1 W fi (R n _ 1 )c(R n _ 1 )f n ^i - 2nR n W if (R n )c i f n - nRlf n /r ex + TcB%_ 1 G nr . 1 f nr -i, (5) 
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Figure 2: The position of exciton condensed phase islands around the laser excitation spot 



R n is the radius of the disk with n excitons, Wfi(R n ) and Wif(R n ) are probabilities for the exciton 
to be captured by the disk and to escape from the disk per unit length of the circle, respectively, 
G n = f G(r)dS n /irRn is the mean value of the exciton pumping over island area, c(R n ) and Cj are 
the concentrations of excitons on the circle of the disk and inside disk = l/s ,s is the area per 
single electron- hole pair in the island). The exciton concentration c(r) differs from n ex (r) presented 
in Fig.l due to presence of islands which perturb the exciton concentration field around them. The 
following condition between transition probabilities W/i and Wif takes place due to the detail balance 
principle 



where W/i(oo) and Wif (oo) are the transition probabilities in the case of straight line boundary 
between condensed and gas phases, Wfi(oo) /Wif (oo) = q/coo, Cqo is the equilibrium concentration 
of excitons for the straight line boundary between condensed and gas phases 




(6) 



oo 



c 10 exp(-<f/KT), 



(7) 



(p is the condensation energy per an exciton, cio = 7(m* kT /2nh 2 , m* is the effective exciton mass,7 
is the degeneracy of the exciton state, a 2 = d[S / kT, a ; is the energy per unit of the disk length. 

The connection between islands occurs due to dependence of the exciton surface concentration 
of considered island c(R n ) in the kinetic equation (j3J) versus presence of other islands. We suggest 
that the distance between disks is larger than the disk radius and the concentration field created 
by some island slowly changes in limits of a size of the considered island. As the result the exciton 
concentration on the disk with m = may by presented in the form 

c ex {R) = n ex {r ) + a K {(R/l) + £ K {2r sin(^/2)//), (8) 

where K (x) is the modified Bessel function, I is the diffusion length of free exciton. Coefficients a M 
are determined by boundary conditions on every disk: the current of excitons to an island should 
be equal difference between number of excitons captured by the island and the number of excitons 
which escape from the island 

27tRD ex ^^- = 2nR(W fl c ex (R) ~ W ifCi ), (9) 

Due to the system symmetry all islands have equal parameters: = ao for every /i. Further we 
introduce the radius distribution function/ (i?) = f n dn/dR = 2irRf n /s . In steady case at n » 1 
the solution of the equation (jlj) has the form 

f(R,N) = f exp(-F(R,N), (10) 

where 



* 2, (nUro) - -r(i-u g+ . „ 



F(R,N) = -4vr / ) ^ ) {-MR, (11) 



o 2v \ n ex {r ) + c 10 e t j-R^i + n G + 

N-l 

a(N) = £ ^o(2f sin(^/2)/[), (12) 

Here we have introduced the dimensionless variables R = Rj yfs^, c(R) = c(R)/s , v = WfiT exy /s^, n G - 

Gs r ex , ip and a are expressed in ([TO]) and (fTTj) in units of temperature. The most probable radius is 
determined from the condition 

dF(R, N) , . 

dfi = °- (13) 
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The mean radius depends on the number of islands is R = R(N). The probability for system to 
have N islands with radiuses Ri,R 2 , ...Rn equals to 

W(N,R u R 2 ,...R N )^exp(-J2F(R t ,N). (14) 

i 

After integrating over the radiuses of the disks we obtain the probability for the system to have N 
islands on the ring 

W(N) = exp(-$(iV)), (15) 

where 

roc _ _ 

$(N) =-N hi z(N),z(N)= exp(-F(R(N),N)dR. (16) 

Jo 

The number of islands is determined by the condition 

The formulae (|13( IT7j) determine the mean radius and number of islands. Since islands of con- 
densed phases derive excitons from the same source two islands can not be situated close one to 
another. Moreover, the distance between islands can not be large because in this case the exciton 
density between them becomes greater than critical value and as the result the new island may ap- 
pear. So, there is specific interaction between condensed phases trough the exciton concentration 
fields. As the result the dependence $(iV) ()17|) has minimum at some value N = N s that deter- 
mines the number of islands on the ring. We have calculated the some parameters of the islands 
on the ring using the equations (|lUH17j) . The exciton parameters and parameters of pumping are 
the same as parameters in Fig.l. The following parameters of the condensed phase were chosen : 
m* = 0.37mo,y? = 10K, a 2 = 18K,s Q = 10~ n cm 2 . It is seen from Fig. 3 that the distance between 
islands (d = 2irr /N s ) increases with rise of temperature, the island radius changes slowly with 
temperature. At considered parameters the condensed phase exists till 2.16K. 

To estimate the island position fluctuation we studied the shift of some island versus the angle tp 
along the ring at fixed position of all other islands (similar to Einstein model for crystal oscillations). 
In this case the function cr(iV) in (|T3*|) depends on ip (cr(N) — > cr(N, if)) and the function f(R, N) — > 
f(R, N, ip) determines the angle distribution function. Using its we calculated the mean fluctuation 
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Figure 3: The dependence of the island radius (R), the distance between islands (d) and the ratio of 
island fluctuation shift to distance between islands (Ad/d) on temperature 

shift of island along the ring, which equals Ad = ro\fp^- It is seen from Fig. 3 that the relative 
shift Ad/d is small (Ad/d << 1 ) and increases only in the vicinity of temperature at which the 
condensed phase disappears. By similar way using (15) the relative fluctuations of the island number 
(AN/N S = \J (N — N s ) 2 /N s ) were calculated. This value is much less than unity. For example, we 
have AN/N S = 7 ■ 1CT 3 at T = 2K. Thus, the studied system has perfect periodicity. These results 
explain the periodical fragmentation observed in |2j. 

We suggested that outside of islands the exciton gas is nongenerated. At considered parameters 
the exciton density at the surface of disks outside of disks equals to 7.2 • 10 8 cm -2 and is significantly 
lower than the critical Bose-Einstain condensation density. Inside of islands the exciton density is 
significantly higher and taking into account the Bose statistic may be important. But presented 
theory did not consider the model of the condensed phase and needs only some parameters of the 
condensed phase: energy and occupied volume per electron-hole pair and surface energy. Moreover, 
the condensed phase may present electron-hole liquid phase. 

In the paper the binding energy per electron-hole pair in condensed phase was chosen (<p = 
10K = 0.8meV). This value gives the critical temperature of fragmentation disappearance (2.1QK) 
close to experimental meaning (~ 2i^).The value of (p is less than width of the fragment emission 
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line (1.3meV [2]) and the emission line spectral positions of free exciton and the fragment coincide 
practically. Even there is very small (much less than width) shift of the fragment emission to blue side 
(See Fig2c of 2\). But mechanisms, that determine the width (scattering by fluctuation potential, 
phonons, between excitons and others), give different contribution to the shapes of emission spectra 
(including to the shift of lines) for the free excitons and the exciton condensed phase. So, if the 
exciton binding energy to condensed phase is less than the bandwidth, the physical information from 
comparison of emission spectra shapes of free excitons and the condensed phase can be obtained only 
together with analysis of line broadening mechanisms for these states. For that the detailed theory 
of the condensed phase is needed. 

Acknowledgment. This research work was supported partially by INTAS grant No. 03-51-5266 
and by the Ukrainian Ministry of Education and Science (project No. 02.07/147). 



Appendix 

To estimate the Van der Waals interaction between two excitons let us use the well known formula 

T r |< 0,0 | U dd | hyk >| 2 , ia , 
Uvdw - - 2^ p ZTp ' ' ' 

where Udd = Pl R f 2 — 3F>1 " Rl pF 2 Rl2 is the operator of the dipole-dipole interaction between excitons, 
Ru is the distance between excitons, 1 , > and \ii,i2 > are the ground and excited states of a 
system with two excitons. The formula Udd can be used if the distance between excitons is larger 
than exciton radius. For calculation we have used the following model. Electron and hole are localized 
in different wells. The distance between wells is equal to S. For estimation we took into account only 
the ground state and two lowest excited states, the dipole transition matrix element of which from 
ground state to excited one does not equal to zero. The wave functions of such type have the form 

|0 >= A exp(— a$r); \la >= A\ a r cos<^exp(— a.\ a r)\ \ \b >= 74ibrsin<^exp(— a^r); (19) 

r is the distance between electron and hole in exciton in the plane of wells, ckq, ax a , are vari- 
ational parameters. The results of calculation of the dipole-dipole and Van der Waals interaction 
(< 0, 0|{7<ta|0, > +U v dw) are presented in Fig. 4 for different values of 5. The energy of bulk exciton 
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Figure 4: The dependence of the sum of the Van der Waals and dipole-dipole interaction energies 
between two excitons in coupled quantum wells on the distance between excitons for the different 
values of the distance between wells 5. E and a are the energy and radius of exciton, respectively, 
in bulk crystal 

is chosen as an unit of energy, radius of bulk exciton (the value of order of 50A°) is the unit of length. 
At the large distance the interaction between excitons is dipole-dipole repulsive. But the attractive 
Van der Waals interaction exceeds the dipole-dipole repulsion at distances less than 3 -j- 6 exciton 
radiuses and total interaction is attractive. The taking into account the other more highest excited 
intermediate states in (jl8p should lead to increasing the presented estimations of Van der Waals 
interaction. So, Van der Waals interaction is really larger. At small distances between excitons the 
exchange interaction becomes important, also the approximation of dipole-dipole interaction becomes 
inapplicable. As the result the total interaction is repulsive one at small distances. 

It should be noted that attractive interaction may be significant at low temperature when con- 
densed phase is created. At higher temperature excitons are distributed uniformly over crystal, 
the mean distance between excitons exceeds the region of attraction and main contribution to the 
exciton-exciton interaction gives the dipole-dipole one. 




9 



References 

[1] A.V. Larionov , V.B. Timofeev. Pis'ma v Zh. Eksp. Teor. 73, 342, (2001). 
[2] L.V. Butov , A.C. Gessard, D.S. Chemla. Nature 418, 751, 2002. 
[3] D. Snoke , S. Denev , Y. Liu , L. Pfeiffer , K. West. Nature 418, 754, (2002). 
[4] L.V. Butov. Sol.St.Comm. 127, 89, (2003). 

[5] D. Snoke ,Y. Liu, S. Denev ,L. Pfeiffer,K. West. Sol.St.Comm. 127, 187,(2003). 

[6] L.V. Butov, L.S. Levitov, A.V. Mintsev, B.D. Simons, A.C. Gossard , D.S. Chemla. Phys Rev. 
Lett. 92, 117404, (2004). 

[7] R. Rapaport, Yang Chen, D. Snoke , Steven H. Simon, Loren Pfeiffer, Ken West, Y. Liu and S. 
Denev. Phys Rev. Lett. 92, 117405, (2004). 

[8] V.I. Sugakov. Fizika Tverdogo Tela 21, 562,(1986). 

[9] V.I.Sugakov. Phase Transition 75, 953, (2002). 

[10] Yu.E. Lozovik , O.L. Berman. Pis'ma v Zh. Eksp. Teor. Fiz 64, 526,(1996). 

[11] V.I. Sugakov. Fizika Tverdogo Tela 46,1455, (2004); Physics of the Solid State 46,1496, (2004). 



10 



